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^-^ — ■ ^ Abstract. Charge-stabilized colloidal suspensions can be conveniently described by 

Q ' formally reducing the macroion-microion mixture to an equivalent one-component 

"^ , system of pseudo-particles. Within this scheme, the utility of a linear response 

approximation for deriving effective interparticle interactions has been demonstrated 

[M. J. Crimson and M. Silbert, Mol. Phys. 74, 397 (1991)]. Here the response 

approach is extended to suspensions of finite-sized macroions and used to derive 

' ^ ' explicit expressions for (1) an effective electrostatic pair interaction between pseudo- 

macroions and (2) an associated volume energy that contributes to the total free 
O I energy. The derivation recovers precisely the form of the DLVO screened- Coulomb 

effective pair interaction for spherical macroions and makes manifest the important 

influence of the volume energy on thermodynamic properties of deionized suspensions. 

^ ' Excluded volume corrections are implicitly incorporated through a natural modification 

of the inverse screening length. By including nonlinear response of counterions to 
macroions, the theory may be generalized to systematically investigate effective many- 
body interactions. 
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1. Introduction 

Charge-stabilized colloidal suspensions, composed of charged niacroions (1 — 1000 nm 
in diameter) and microions (counterions and salt ions) suspended by Brownian motion 
in a molecular fluid, occur in a variety of natural and engineered forms |1]]. Common 
examples include clay minerals, paints, inks, and detergents (micellar solutions), in 
which repulsive electrostatic interactions promote stability against coagulation induced 
by van der Waals attractive forces. Somewhat more exotic are suspensions of synthetic 
latex or silica spheres whose near monodispersity facilitates self-assembly into crystalline 
lattices 0. Aside from serving as well-characterized models for fundamental study, 
synthetic colloids exhibit unique optical properties that underly several emerging 
technological applications, such as nanosecond optical switches [^ and photonic band 
gap materials Q. 

Accurate prediction of the physical properties of colloidal matter relies on a 
fundamental understanding of interparticle interactions. The first quantitative account 
of electrostatic interactions was achieved by Derjaguin, Landau, Verwey, and Overbeek 
(DLVO) 0. Based on the Poisson-Boltzmann equation for the electrostatic potential |l[], 
the DLVO theory portrays the bare Coulomb interactions between macroions as 
effectively screened by a surrounding atmosphere of microions. The resulting screened- 
Coulomb pair potential has been a valuable cornerstone of colloid science for half 
a century. Nevertheless, experimental evidence for apparent long-range attractions 
between macroions Q has contributed to renewed interest in colloidal interparticle 
interactions. 

An explicit description of the multi-component mixture of macroions, counterions, 
salt ions, and solvent molecules clearly poses a formidable challenge. Consequently, 
interactions in such complex systems are usually treated at the level of effective 
interactions. Tracing out from the partition function statistical degrees of freedom of 
all but a single component reduces the problem to that of an equivalent one-component 
system of "pseudo-particles" governed by an effective thermodynamic-state-dependent 
interaction [0. A range of theoretical and numerical methods have been developed 
to investigate effective interactions in charge-stabilized colloids. Strategies deployed 
to date include Poisson-Boltzmann cell models P], Monte Carlo simulation |^, p!0[] , 
ab initio simulation [11|, and density-functional theory [12, 13, n, 15]. Recently, 



Silbert and coworkers |T6|, [T^ proposed an approach motivated by analogies between 



charged colloids and metals. With the correspondences (counterion <-;> electron) 
and (macroion •«-»• metallic ion), the procedure is closely akin to the pseudopotential 



theory of metals |]T8|, which can successfully account for thermodynamic properties 
of simple metals [jTU[. Performing a classical trace over microion degrees of freedom 
and describing the electrostatic response of the microions to the macroions within 
second-order perturbation theory, results in an effective interaction between pseudo- 
macroions and an associated volume energy that contributes to the total free energy. 
The importance of including the volume energy in calculating thermodynamic properties 
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of charge-stabilized colloidal suspensions has been emphasized recently by a number of 
authors [0, 0, |T|, |1|, ^|. 



The main purpose of this paper is to develop further the response approach to 
charge-stabilized colloids, extending it, in particular, to suspensions of finite-sized 
macroions. The proposed extensions enforce exclusion of microions from the macroion 
hard cores and explicitly take into account the volume excluded by the macroions to 
the microions. A second goal of the paper is to establish a framework for generalizing 
the theory to include nonlinear response to allow investigation of effective many-body 
interactions. In the next section, following a brief review of the response approach, 
the extensions are outlined. Section 3 presents the main results - obtained within a 
linear response approximation - for an effective pair potential acting between pseudo- 
macroions and an associated volume energy, both of which consistently incorporate 
excluded volume effects. Finally in Sec. 4, we discuss implications of the results for 
thermodynamic properties of charge-stabilized colloids and prospects for generalizing 
the theory beyond linear response. 

2. Theory 

2.1. The Model 

The theory described below is based on the "primitive" model, wherein for simplicity 
the solvent is treated as a uniform dielectric continuum. To simplify notation, it 
is furthermore assumed that counterions are the only microions present (deionized 
suspension). The general case of finite salt concentration will be addressed 
elsewhere 0] . The model system then consists of iVm charged hard-sphere macroions of 
diameter a and charge —Ze (e being the elementary charge) and A'^c point counterions 
of charge ze suspended in a uniform fluid medium characterized entirely by a dielectric 
constant e. Each macroion is assumed to carry a fixed charge, uniformly distributed 
over its surface. Charge fluctuations are thus implicitly ignored in the model. Described 
statistically by a canonical ensemble, the system occupies a fixed total volume V at 
temperature T. For a given number of macroions, global charge neutrality constrains 
the number of counterions by the condition zN^ = ZN^. 

Denoting macroion and counterion positions by Rj and Tj, respectively, the full 
Hamiltonian of the system can be expressed in the form 

H = Hja + Hc + -f^mc, (1) 

with the individual terms to be specified below. The first term on the right side of 
Eq. (|l]) is the bare macroion Hamiltonian, given by 

1 ^"^ r 1 

H^ = K^ + -Y^ [vns{\ni - Kj\) + t;„,„,(|R. - R,-|)J, (2) 

i,j = l 

where K^ is the kinetic energy of the macroions, fHs(|Rj ~ R-jl) is a hard-sphere pair 
interaction between the macroion cores, and fmm(^) = Z'^e'^/er is the bare Coulomb 
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interaction between a pair of macroions whose centres are separated by a distance r > a. 
The second term in Eq. (|1|) is a counterion Hamiltonian, taking the form 

He = Kc + ^E^cc(|r.-r,|) + $:^t;Hs(|r.-R,|), (3) 

^ ij = l i = lj = l 

where K^ is the counterion kinetic energy, fcc(^) = z^e^ /er is the Coulomb interaction 
between a pair of counterions, and fHsdrj ~ R-il) is the hard-sphere interaction between 
a point counterion and a macroion core. Finally, the third term in Eq. (0) is the 
electrostatic interaction energy between the macroions and counterions: 

^mc = EE^mc(|r.-R,|), (4) 

where fmc(^) denotes the macroion-counterion electrostatic pair interaction. Outside the 
macroion core radius, fmc(^) has the Coulomb form. Inside the core, however, v-mdj") is 
not uniquely defined. Thus, following van Roij and Hansen |]T^, we choose v^cij") to be 
a constant for r < a and take 

^2 



^#^, r > a/2 

-a, r < (t/2 



ea/2 

where the parameter a will be specified (Sec. 3.1) to ensure that the counterion density 
vanishes within the core, as physically it must. 

2.2. Reduction to an Equivalent One-Component System 

Having defined the Hamiltonian, we now turn to a statistical mechanical description, 
with the ultimate aim of calculating the free energy of the system. The partition function 
is given by 

Z^ = {{exp{-H/k^T))X, (6) 

the angular brackets symbolizing classical traces over macroion or counterion degrees 
of freedom. Following standard treatments developed first in the context of simple 
metals |IU| , ^ , we proceed by reducing the two-component mixture of macroions and 



counterions to an equivalent one-component system. The reduction is achieved by 
performing a restricted trace over counterion coordinates, while keeping the macroion 
coordinates fixed. Thus, without approximation in this purely classical system, 

Zm = (exp(-i7eff/fcBT))^, (7) 

where H^g = H^ + Fc is the effective Hamiltonian of a one-component system of pseudo- 
macroions, and where 



F, = -kBT\n[exp[-{H, + H^,)/kBT\)^ (8) 

may be physically interpreted as the free energy of a nonuniform gas of counterions in 
the presence of macroions fixed at positions Rj. In general, the counterion free energy 
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is a complicated many-body function of the macroion positions. Progress can be made, 
however, by formally adding to and substracting from H a term, E]^, representing the 
energy of a uniform compensating negative background. Then Ec may be expressed in 
the form 



F, = -A;BTln(exp[-(//^ + if;j/A;BTj)^, (9) 

where H'^ = Hc + E^, and H^^^ = H^ac — E^,. The advantage of this manipulation is that 
H'^ is simply the Hamiltonian of a classical one-component plasma (OCP) of counterions 
in the presence of neutral hard-sphere macroions. 

Now the counterions, being excluded by the hard cores of the macroions, occupy 
the free volume Vf = V — N^{n/Q)a^, i.e., the volume not occupied by the macroion 
cores. The average effective density of counterions is therefore given by ric = Nc/Vf = 
n^°y{l — r]), where n^°^ = Nc/V is the nominal counterion density and r] = {V — Vi)/V 
is the macroion volume fraction. An important question now arises. What volume 
should the background occupy? In order that H'^ truly be the Hamiltonian of an 
OCP, the background and counterions clearly must occupy the same volume. In fact, 
were the background to occupy a different volume {e.g., the total volume V) then 
the effective Hamiltonian would contain terms that are formally infinite (see below), 
although identically cancelling, associated with the long-range Coulomb interaction. 
Thus, the background is taken to be excluded - along with the counterions - from the 
macroion cores, its density equaling the effective counterion density ric- 

The background energy is then given explicitly by p2| 



where 



^b = - -nl / dr / dr'^^ = - -N,n,v,M. (10) 

2 Jvt Jvt e|r — r I 2 

«cc(0) = / dr^ = lim(^) (11) 

Jvi er k^Q\ ek ' 

is the A; — i> limit of the Fourier transform of fcc(^)- The infinity arising from E\^ will 
be seen below to be formally cancelled by an identical infinity in H^c- 

2.3. Linear Response Approximation 

Thus far, the theory is exact, within the primitive model. The challenge lies ahead in 
calculating the counterion free energy [Eq. (||)]. One proposed approach [J12l invokes 



density-functional theory to approximate E^., regarded as a functional of the counterion 
density, by expanding in a functional Taylor series about a uniform counterion OCP. 
An alternative strategy [|^, 0, inspired by the pseudopotential theory of metals, is 



to formally regard H'^^ as an "external" potential acting upon a counterion OCP and 
approximate E^ by perturbation theory. Following the latter strategy [p^ , we write 



Fc 



Focp + f'dX{H'^,),, (12) 

"^ U 
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where 

i^ocp = -kBTlniexpi-H'Jk^T))^ (13) 

is the free energy of the reference counterion OCP, occupying a volume Vf, in the presence 
of neutral hard-sphere macroions. The integral over A in Eq. (0) physically corresponds 
to an adiabatic charging of the macroions from neutral to fully-charged spheres. The 
ensemble average ( )^^ represents an average with respect to the distribution function of 
a system whose macroions are identically configured but carry a charge XZ. 

Further progress is facilitated by expressing (H^^)^ in terms of Fourier components 
of the macroion and counterion densities and of the macroion-counterion interaction, 
according to 

(^mc)A = 77 E ^mc(fc) (pc(k));, p„,(-k) + - limfe^c(fc) {pc{k)) Jm{-k)] - Ey,. (14) 

Evidently (H^^)^ depends through Pc(k) upon the response of the counterions to the 
macroion charge density. Regarding the macroion charge as imposing an external 



potential on the counterions, the counterion density may be expressed in the form p3 

Pc(k) = X^^H^)^mc(fc)Pm(k) + — X!x^^Hq>k-q){)mc(g)^)mc(|k-q|)pm(q)Pm(k-q) 

>/f q 

+■■■, (15) 

where x is the i*^ member of a hierarchy of response functions of the reference 
counterion OCP. Here, as in ref. [1^, we adopt the simplest nontrivial approximation 
and assume that the counterions respond linearly to the macroion charges. Although 
its range of validity is uncertain, linearization is expected to be justified for sufficiently 
dilute suspensions and weakly charged macroions. Thus we take 

(Pe(k));, = x{k)XVrnc{k)prn{k), k ^ 0, (16) 

where x{k) = X^^Kk) is the linear response function. Note that for /c = there is no 
response, since Pc(0) = A^c is fixed by the number of counterions. Substituting Eqs. (^) 
and (|T^) into Eq. (p!2| ) and integrating over A, the counterion free energy is given to 
second order in the macroion-counterion interaction by 



1 



-^c = -^ocp + Tprr J2 X{k) [vmc{k)] pm(k)pm(-k)+nclim N^v^dk) + ^Vcc{k) .(17) 



A^. 



Correspondingly, the effective Hamiltonian takes the form 
1 ^™ 1 

HeS = A'm + - E ^Hs(|H.i " R-jl) + 7^77 E ^mm(/i;) Pm(k)pm(-k) - N^ 

+ ^ocp + ttttExI^) [Omc(A;)]^Pm(k)pm(-k) 

+ Tie lim — 7r^x{k) [Vrac{k)f + N^Vmc{k) + -^Vcc{k) 



2 
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Notice, however, that Eq. ( ]T8| ) may be restructured and written in the form 

^ «., = ! ^^f k 



1 JVm 

= ^m + ^ E [^HS(|R^ - RjI) + t^eff(|Ri - Ril)] + ^o, (19) 

where 






Vesik) = V^mik) + Vindik) (20) 

may be interpreted as an effective electrostatic pair potential between pseudo-macroions, 
being the sum of the bare Coulomb potential and an induced potential 

i)ind(A;) = x{k)[Vrnc{k)]\ (21) 

The final term in Eq. (0), 



N„ 



z . ,,. . .,. Z 



Eo = Focp + -— limt;ind(r) + ricN^n lim -j-TzVmd{k) + Vr^cik) + -—Vcc{k) , (22) 

Z 1 — >0 fc-^OL ZZ Zz J 

is the volume energy, which is a natural and inevitable consequence of the reduction 
to an equivalent one-component system. Although having no explicit dependence on 
the macroion positions (see below), Eo evidently depends on the average density of 
macroions and thus can make a significant contribution to the total free energy of the 
system. It must be emphasized that the above expressions for the effective pair potential 
and the volume energy are identical to expressions derived from the pseudopotential 



theory of metals |T9|, 0, ^ if one substitutes for Fqcp and xik), respectively, the 
energy and linear response function of the uniform electron gas, and for v^c{k) the 
electron-ion pseudopotential. 

To summarize thus far, starting from the primitive model of charge-stabilized 
colloids, formally reducing the two-component macroion-counterion mixture to an 
equivalent one-component system of pseudo-macroions, and applying a linear response 
approximation to the counterion density, we have obtained expressions for both an 
effective electrostatic pair interaction [Eqs. (|20| ) and (0)] and an associated volume 
energy [Eq. (|22|) ]. Practical calculations still require specification of (1) the reference 
OCP free energy Fqcp, (2) the OCP linear response function x{k), and (3) the macroion- 
counterion interaction 'Omc(^)- Below each of these is considered in turn. 

It is important first to note that by associating the hard-sphere part of the total 
macroion-counterion interaction with the counterion Hamiltonian [Eq. (|^)] - necessary, 
since response theory does not apply to hard-sphere interactions - the reference OCP 
is confined to the free volume between the macroion cores. As a consequence, the 
OCP is not strictly uniform since, in principle, the boundary conditions may induce 
nonuniformity. Determining the free energy of such a system in general poses a nontrivial 
task. In practice, however, counterion densities are usually low enough that the OCP 
may be assumed to be essentially uniform, except perhaps near contact with a macroion 
surface. 
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Now, for typical macroion charges and concentrations, the OCP is so weakly coupled 
(unlike its electronic counterpart in metals) that its free energy is dominated by the 
ideal-gas entropic component. Therefore, ignoring correlations between counterions [0, 
13, E0|, an accurate approximation is 



Focp ^ kBTNMn^A')-!] = fcBTiVjln(^^(^^^^) - ll , (23) 

L J L \ 1 _ 7^ / J 

where A is the thermal de Broglie wavelength and rim = N^/V is the average number 
density of macroions, the last equality following from the constraint of global charge 
neutrality. 

The linear response function is directly related to the corresponding static structure 
factor S{k) via 

X{k) = -f3n^S{k) = - ^\ (24) 

where /3 = l/k^T and c{k) is the Fourier transform of the direct correlation function 
c{r). Specifying x{k) is therefore equivalent to specifying c{k). For a weakly coupled 
OCP, a convenient and reasonable approximation for c{r) is given by the mean spherical 
approximation (MSA). This amounts to setting c{r) equal to its asymptotic (r -^ oo) 
limit c{r) ~ —l3vcc{r) for all r. As a result, 

c{k) ^ -/?f)cc(fc) = -^^- (25) 



Substitution of Eq. (H) into Eq. (H) then yields [|T6 

1 + K'^/k^ 



Xik) = - T^^, (26) 






where 

efceT ^ ~ ^l-r])ekBT^ 
As will be seen below, the parameter k, plays the role of an inverse screening length in 
the counterion density profile and in the effective pair interaction. 

Finally, specifying the macroion-counterion interaction amounts to determining the 
value of the parameter a in Eq. (^) that will ensure a vanishing counterion density inside 
the macroion cores. This in turn requires a calculation of the real-space counterion 
density profile, the details of which are described in the next section. 

3. Results 

3.1. Counterion Density Profile 

The real-space counterion density profile Pc(r) may be determined from Eqs. (^, ([T6|), 
and (|26|) . First, Fourier transforming Eq. (|^) yields 

AnZze'^r, , ,, , ^ sin(A;cr/2)i , , 

Vmc{k) = - ^^^ [il-a)cos{ka/2) + a ^^ ' ^ ]. (28) 



Effective Interactions in Colloidal Suspensions 

Next, substituting Eqs. (|2^) and ( pH]) into Eq. (|1^) gives the intermediate result 



Z / K 



2 



sin(A;cr/2)- 



Pc(k) = j{j^) [(!-«) oos{ka/2) + a lly^'^ ] Eexp(.k.R), (29) 

where the sum is over the positions R of the macroions. For simphcity, we consider 
the density profile around a single macroion located at the origin (R = 0), assuming 
all other macroions to be far away [nR ^ 1). This is equivalent to retaining only the 
R = term in the summation, in which case Pc(k) is a function only of k and Pc(r) is a 
function only of the radial distance r. Now inverse Fourier transforming Eq. (E9|) yields 




X , . ,„x , smh(Kcr/2)i exp(-Kr) ^ ,„ 

a)cosh(Ka/2) + a ^ ,' ' — ^-\. -, r>a2 

Pcl'^J - <; V .2 , . _ sinh(fi:r) ,^ ^'^^> 

Vanishing of Pcir) for r < o j2 is evidently ensured by setting 

KO 12 , , 

a = '-—. 31 

1 + kct/2 ^ ^ 

Finally, substituting this expression for a back into Eq. (^) gives the result 

Pc(r) = -f -^— -^, r>a/2, (32) 

which is automatically normalized to the correct number of counterions per macroion 
[Zjz). This expression for the counterion density profile around a single macroion 
is recognized to be of precisely the same form as the Debye-Hiickel expression for 
the density of electrolyte ions around a macroion [Q], where k is the inverse Debye 
screening length. A notable distinction lies, however, in the definition of k. Whereas 
our K [Eq. (^)] depends on the average effective counterion density ric in the volume 
unoccupied by macroions, the Debye-Hiickel k depends rather on the nominal bulk 
density of electrolyte ions. The importance of redefining the usual k in this way - a 
result that emerges naturally from the response approach - has been stressed also by 
Russel and coworkers p5| . 

In passing, two remarks are in order. First, determining the counterion density 
profile in the presence of two or more closely spaced macroions will evidently require 
a more general form for the macroion-counterion core interaction than the simple 
constant chosen in Eq. (|^). Second, it may be instructive to compare the macroion- 
counterion interaction with its metallic counterpart, the electron-ion pseudopotential. 
A popular and successful form of the latter is the empty- hole pseudopotential [^ , which 
is Coulombic at long range but precisely zero inside a certain core radius. In contrast, 
setting fmc(^) to zero for r < oj2 would result in a nonvanishing counterion density 
inside the macroion cores. A distinction between the metallic and colloidal cases lies 
in the fact that, while counterions are strictly excluded from macroion cores, electrons 
may at least partially penetrate metallic ion cores. 
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3.2. Effective Pair Interaction and Volume Energy 



10 



We are now in a position to derive the main results of the paper. Considering first 
the effective electrostatic pair interaction between pseudo-niacroions, we proceed by 
substituting Eq. (|3lD into Eq. (p8l), obtaining 



VracK,!^) 



A-KZze^ 



(r 



f 



cos(A;cr/2) + K, 



sin{ka/2) 



Next substituting Eqs. (pBD and 
. ,,, 27rZV/ 1 ^2 

Vmd{k) = - 



+ Ka/2 



into Eq. (|21| ) yields 



k 



eP Vl + Kcj/2/ VP + k2 



,, , „ sm(ka) ol — cosf/ccr) 
l+cos{ka)+2K — \ — -+K^ ^ ' 



k 



fc2 



(33) 



.(34) 



Fourier transformation of Eq. 
result 



is a straightforward, if tedious, calculation, with the 



^2^2. gW2 ^2 g-^r- ^2^2 



^'md(r) 



(r 



e Vl + K^ 



er ' 



"2iF 



.f L 



(2 + Ka)Kr — ^K^r'^ 



r > a 
r < a. 



(35) 



+ Ka/2j 

Finally, substituting Eq. (^) into the Fourier transform of Eq. (^), we obtain an explicit 
expression for the real-space form of the effective electrostatic pair potential: 

^2g2 ^ ^Ku/2 



Vesir) = v^^{r) + Vindir) 



2 e" 



r > a. 



(36) 



e ^1 + Ka/2^ r 

This result is identical in form to the electrostatic part of the familiar DLVO effective 
pair potential 0, which is usually derived from the Poisson-Boltzmann equation. The 
only difference between our potential and the DLVO potential lies in the definition of k 
[Eq. (p?!)], which here involves the effective counterion density n^. 

The volume energy now may be explicitly determined from Eq. (|22D. It follows 
immediately from Eq. (|35D that 

K 



2^2 



liniUind(r) 

1 — >0 



Z'e 



from Eq. 



limT)ind(/i;) 



that 



Zn2^ , ^ 47rZV 

- ^cc(0) + 



Ka/2'' 



nZ^e^a'' 



en" 



1 + Ka/2 



and from Eq. (|33D that 



\imv^c{k) 



Z ^ ,^, -nZze^ 
z 2e 






^ 3 12 ^ 



:i + -«:a). 



(37) 



(38) 



(39) 



Substituting Eqs. (|3 
the volume energy: 

E. -- 



+ Ka/2' 
and (|39| ) into Eq. (p2D, we obtain the following result for 



OCP 



iV„ 



Z'e 



2^2 



K 



A^r 



Z^bT 



(40) 



2e l + Ka/2 '"' 2^ 

The first term on the right side of Eq. (^OD is the OCP free energy, discussed in 
the previous section. The second term, which depends implicitly on the macroion 
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density through the parameter k, may be given a physical interpretation as one half 
the electrostatic energy associated with a single pseudo-macroion, composed of a 



macroion surrounded by its own screening cloud of counterions [12, Q. The final 
term, corresponding to the k ^ limit in Eq. (^2]), contributes a density- independent 
constant to the free energy per macroion and hence has no influence on thermodynamic 
phase transitions at zero salt concentration. At finite salt concentration, however, the 



corresponding term cannot be ignored |I3|, |1J, ^ 

4. Discussion and Conclusions 

It is important to point out some limitations of the theory and the results presented 
above. First of all, the assumption of linear response of the counterions is strictly valid 
only for dilute suspensions of weakly charged macroions. Whether the linear response 
approximation remains valid at higher concentrations - in particular, concentrations 
for which excluded volume effects begin to play a role - is an interesting and open 
question. In this regard, it may be worth noting that Poisson-Boltzmann cell model 



calculations M and ab initio simulations |11| do support the general form of the screened- 



Coulomb pair potential at appreciable concentrations, albeit with renormalized DLVO 
parameters. Secondly, by considering only electrostatic and hard steric interactions, 
and ignoring short-range interactions between counterions and macroion surfaces, the 
theory cannot address the possibility of condensation of counterions onto the macroions 
and the consequences for effective macroion charges. Finally, no account is taken of 
correlations between charge fluctuations, either on the macroion surfaces or in the 
density distributions of counterions surrounding neighbouring macroions. For spherical 
macroions, Monte Carlo simulations and cell model calculations |27| suggest that such 



correlations make only a small contribution to the total free energy, at least at low salt 
concentrations. On the other hand, for rod-like macroions, correlated charge fluctuations 
may play a more significant role 1 2^ . 



In summary, by applying a linear response approximation to the counterions 
surrounding charged monodisperse hard-sphere macroions in a colloidal suspension, we 
have derived two main results, namely an effective electrostatic pair interaction Ves{f) 



[Eq. (^)] and an associated volume energy Eq [Eq. (^Ol)]. The total free energy of 
the system is ultimately the sum of the volume energy - whose physical origins are 
the counterion entropy and the macroion-counterion interaction energy - and the free 
energy of the equivalent one-component system of pseudo-macroions interacting via 
their effective pair interaction. Our expression for Vcs{r) confirms that a linear response 
approximation, combined with the MSA for the response function, yields the familiar 
DLVO form of pair potential for spherical macroions, prefactors included [^Sj. This is 



not surprising, given that the DLVO potential also may be derived by linearizing the 
Poisson-Boltzmann equation. At the same time, however, our derivation indicates how 
excluded volume corrections may be incorporated through the density dependence of the 
inverse screening length k, by substituting for the nominal counterion density n'f^ the 
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effective density Uc of counterions occupying the free volume between macroion cores. 

Our expression for the volume energy, which quantitatively exhibits the dependence 
on the macroion density, confirms the necessity of including Eo in calculating 
thermodynamic properties from the free energy [|1^]. In particular, the phase behaviour 
of deionized suspensions of highly charged macroions have been shown to depend 
sensitively on the volume energy [T^, ^, H, 2D|. Bulk pressure and elastic constants are 



expected to be similarly sensitive. Furthermore, our expression for E^ is consistent with 
that obtained by van Roij et al [jTsf from an alternative density-functional approach. The 



sole distinction is that our expression, which involves the effective counterion density ric, 
incorporates excluded volume corrections, at least in an approximate fashion. Although 
not likely of significant consequence at the small volume fractions considered in Refs. 
and |jl3|, such effects may become important at higher concentrations [pOf . 



In conclusion, the linear response approach of Silbert et al [0, |13 offers a powerful 
tool for investigating effective electrostatic interactions in charge-stabilized colloidal 
suspensions. As demonstrated here, a consistent extension to finite-sized macroions 
leads directly to (1) an effective pair interaction between pseudo-macroions having 
precisely the DLVO screened- Coulomb form, but with a modified inverse screening 
length that incorporates excluded volume corrections, and (2) a density-dependent 
volume energy that can make a significant contribution to the total free energy of 
salt-free suspensions. By including the next higher-order response function [x^^^^ in 
Eq. (0)], the approach can be straightforwardly generalized to include nonlinear 
response of microions and thereby used to assess the importance of effective three- 
body interactions |]TT|. This is equivalent to approximating the counterion free energy 
[Eq. (|12|)] in perturbation theory to third order in the macroion-counterion interaction. 
An expression for an effective triplet interaction already has been derived from a 



density-functional approach |30|. It remains, however, to analyse the corrections that 
are entailed both to the volume energy and to the effective pair interaction, and to 
explore the implications for thermodynamic properties. Experience from ah initio 
simulations [0 and from the realm of metals [0] suggests that many-body effects 
become significant at sufficiently high densities. Further outstanding issues are whether 
in bulk the effective pair interaction always retains its screened- Coulomb form, and 



whether near a boundary the interaction can ever become attractive |15|. It is hoped 



that in future the response approach may help to resolve these important issues. 
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